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The three dimensional harmonic oscillator model including 
a cranking term is used for an energy variational calculation. 
Energy minima are found under variation of the three oscilla- 
tor frequencies determining the shape of the system for given 
values of the three components of the rotational vector which 
determines the orientation of the angular momentum in the 
intrinsic frame. Tilted rotations are established by numer- 
ical means for triaxial nuclei. The onset of tilted rotations 
is related to an instability of the mean field as verified by 
the random phase approximation. Analytic expressions are 
derived for the critical rotational frequencies associated with 
the mean field instability. 



PACS numbers: 21.60.Ev, 21.60.Jz, 21.10.Re, 03.75.Fi 



I. INTRODUCTION 

Numerous experimental observations of well developed 
rotational bands implying A/ — 1 sequences (see for a 
review Q) raise the question about mechanisms creat- 
ing such collective states. The conventional description 
within the Cranking Shell Model (CSM) (|] suggests the 
A/ = 2 sequence of rotational states. As is well known, in 
the CSM the axis of rotation coincides with the principal 
axis of the density distribution of the nucleus (principal 
axis cranking or PAC). As a result, the CSM Hamilto- 
nian adheres to the D2 spatial symmetry with respect to 
rotation by the angle tt around the rotational axis, for 
instance, the x— axis. Consequently, all rotational states 
can be classified by the quantum number called signature 
r — exp(— iTTQ;) leading to selection rules for the total an- 
gular momentum I = a + 2n, n = 0, ±1, ±2 .... In par- 



ticular, in even-even nuclei the lowest rotational (yrast) 
band characterised by the positive signature quantum 
number r = 1 (a = 0) consists of even spins only. 

It was suggested that non-principal axis rotations 
which invalidate the concept of signature symmetry is 
responsible for rotational bands with A/ ~ 1 sequences. 
This idea has led to considerable experimental and theo- 
retical efforts during the past few years . Related ques- 
tions have been addressed recently in connection with 
effects of superfluidity in a trapped Bose-Einstein con- 
densate (EEC) g. 

For a nucleus, the non-principal axis rotation is con- 
sidered stationary. For the BEC, a precession of the 
symmetry axis is assumed for the condensate around the 
symmetry axis of the confining trap with the precession 
frequency being fixed by the angular momentum carried 
by the vortex line. A quantised vortex is aligned along 
the symmetry axis of the condensate carrying one unit 
h of angular momentum per particle. Similarly, the pre- 
cessional motion of the principal axes in the nucleus cre- 
ates vibrational (wobbling) excitations above the yrast 
line 1^-0]. Notice that in this case the wobbling excita- 
tions arc quantum mechanical fluctuations of the angu- 
lar momentum around a fixed principal axis of the mean 
field generated by the CSM |^,^. These excitation carry 
one unit of angular momentum, similar to the vortex in 
the BEC, and corresponding states may be associated 
with odd spin partners of the gamma band in a rotat- 
ing nucleus at low angular momenta It has been 
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demonstrated |9| that at high rotational frequency the 
wobbhng excitations become particularly low. Similar 
to the transition from spherical to deformed nuclei |^] 
when quadrupole excitation energies tend toward zero, 
an analogous mechanism is expected to prevail in that for 
vanishing wobbling excitations the nucleus undergoes a 
transition to tilted rotation as the energetically favoured 
state. The analysis of the link between tilted rotation 
and wobbling motion is one major aspect of the present 
paper. 

In fact, a planar tilted rotation effecting the transi- 
tion from the superfluid to the normal phase has been 
discussed within the phenomenological two-phase model 
The first microscopic attempt to describe a non- 
principal axis rotation has been carried out in using 
the triaxial cranked oscillator model (TCO) with con- 
stant angular frequency (see for a review of latest studies 
101). However, the TCO calculations employ a fixed tri- 
axial deformation. The violation of the self-consistency 
between rotation and evolution of shape parameters is 
a serious drawback of this model. In the present pa- 
per we investigate the occurrence of tilted rotation using 
the self-consistent triaxial cranked oscillator model. The 
model guarantees full self-consistency at all rotational 
frequencies and all particle numbers. A major outcome 
of the present paper is the result that, in even-even nu- 
clei, tilted rotations occur if and only if the nucleus has 
a triaxial shape in its ground state. The pairing inter- 
action could have an effect upon the conclusions drawn 
in the present paper. However, its role is less important 
at high spins where tilted rotations should be observed 
0. On a technical footing, it involves additional vari- 
ational parameters which complicates the minimisation 



procedure. Therefore, we exclude in our model the pair- 
ing interaction. The transparency of the model which 
provides a general tendency of the phenomena and the 
succinct results compensate for this shortcoming. 

In Sec. II we review the main features of our model. 
The numerical results are presented in Sec. HI. The re- 
sults are discussed in Sec. IV and a short summary in 
Sec.V. Technical details are deferred to Appendices. First 
preliminary results have been reported in ||l^ , [l3|] . 

II. THE MODEL 

The many-body Hamiltonian (Routhian) in the rotat- 
ing frame is given by 

N 

H ^^(hoii) -n-i{i)) ^ Ho-n- L (1) 

where the single particle triaxial harmonic oscillator 
Hamiltonian Hq is aligned along its principal axes and 
reads 




The rotational vector of the cranking term has the com- 
ponents (ila;, rij,, fi^) = {sin 6 cos (j), sin 6 sin (p, cos 9). 

In Appendix A it is shown how this Hamiltonian is cast 
into the form 

N 3 

^ = + 1/2), (3) 

where the normal mode energies are obtained from 
a cubic equation and the normal mode operators Qk are 
linear superpositions of the pi and Xi . In this way we have 
the exact single-particle energies and wave functions at 
our disposal. Note that, for fl ^ and in particular 
for a tilted rotation, the normal mode operators Qk are 
superpositions of all three Cartesian components of the 
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Pi and Xi. Furthermore, the eigenmode energies Ek are 
real only for a certain range of f2 < rJmax (see Appendix 
B). We will restrict ourselves to this physical range of 
values of the rotational frequency. 

The knowledge of the eigenmodes and eigenfunctions 
enables us to calculate energy contours of the total energy 



EiYli + E2T,2 + £-3^3 



(4) 



in, say, the lu^ — ujy plane with the particle number N 
and the rotational vector Q as parameters. The configu- 
rations are determined by = (rife + 1/2) j where the 
occupation numbers are the eigenvalues of Q^Q^, they 
take the values 0,1,2,.... For the PAC rotation around 
the x— axis the notation Si = T^x, S2 = 2+, E3 = E_ is 
used jl^ (we chose ujx > t^y > i^z, see also Appendix B). 

The ground state is determined by filling the single- 
particle levels from the bottom. We take care of the 
particle spin only in obeying the Pauli principle which 
allows two particles in one level; also we consider only 
one kind of nucleons, protons or neutrons. It is clear 
that different sets of normal modes yield different sets of 
occupation numbers. A pronounced shell structure can 
exist only for special sets of normal modes. Correspond- 
ingly, all spherical shells will be strongly mixed in our 
model, in particular the AiV = 2 mixing is taken into ac- 
count exactly. In this way we go beyond the calculation 
of |,|l| where mixing between major oscillator shells is 
ignored. 

The normal modes depend on the three components of 
the rotational vector and on the harmonic oscillator fre- 
quencies. From our assumption that the system adjusts 
itself under the influence of the rotation by minimising 
-Etot , a change of the magnitude or direction of the rota- 
tion leads to a corresponding change of the effective mean 



field potential which is given by the oscillator frequencies. 
In other words, for a given rotational frequency, we must 
seek the minimum of i?tot under variation of the oscillator 
frequencies. The variation cannot be unrestricted as the 
confining potential encloses a fixed number of particles, 
and assuming that the particle density does not change 
we are led to a fixed volume constraint 



UxLOyUj^ = ujo^ = {Alf/N. 



(5) 



The volume constant is chosen in units where frequencies 
are given in MeV and lengths in Fermi. Introducing the 
Lagrange multiplier A, we solve the variational problem 



S{{g\H\g) - XuJxUJyUJz) = 



(6) 



where \g) denotes the ground state as described above. 

From Eq.(^) we obtain, after differentiation with re- 
spect to the oscillator frequencies and using Feynman's 
theorem \M 



^{9\H\9) = {9\^\9), 
diOk diOk 



the self-consistency condition 



u:l{g\x^\g)^ul{g\y^\g)^u^l{g\z^\g) 



(7) 



(8) 



which must be obeyed at the minimum of E'tot- Notice, 
that Eq.(^) becomes a general result which includes as 
a particular case the cranking harmonic oscillator with 
the PAC rotation around the x-axis |l^. The condition 
serves as an acid numerical test. Note that it is fulfilled 
irrespective of the minimum being local or global. 

For f2 = we recover the well known magic numbers of 
the triaxial harmonic oscillator for one type of particles 
(proton or neutron) D at = 20, 40, 70, 112, . . . which 
manifest themselves in our procedure by the spherical 
symmetric solutions ujx = ujy = lOz — wq. Between these 



magic numbers the well known prolate or triaxial near 
prolate and oblate or triaxial near oblate solutions arc 
found, where a (near) prolate solution (uj^ — LOy > t^z) is 
obtained for a shell which is less than half completed and 
a (near) oblate solution iuj^ — < i-^z) foi' a more than 
half completed shell. In Fig.l the values of (x^), (y^) 
and {z'^) are plotted versus N. From the figure we see, 
that, for = 0, some shapes are axially symmetric while 
others are not. Switching on the rotation Q ^ a, similar 
picture is obtained as long as the rotational frequency 
is kept at a reasonably low value (fi < 1.0 for iV < 24 
and n < 0.3 for N > 40). The corresponding minima 
of the total Routhian are a function of il, i.e. of fi, 6 
and (j). Our interest is focused upon the lowest minimum 
under variation of 9 and (p for fixed fi, i.e. dEtot/dO — 
dEtot/dip — 0. If such minimum is found for 7^ we 
refer to this solution as to a tilted solution. 

III. NUMERICAL RESULTS 

As is well known ||l^, the PAC rotation of prolate 
systems leads to triaxial shapes. With increasing ro- 
tational speed the change of the shape leads eventu- 
ally, for a critical frequency ^l^jt, to an oblate shape 
with the rotational axis coinciding with the symmetry 
axis. Using the minimising procedure of our model con- 
firms this fact in principle. However, our model allows 
an easy assessment as to whether the minimum is lo- 
cal or global. It turns out that, for N < 56, the min- 
imum found at ^l^J^^ is no longer a global minimum, 
since somewhere in the interval [0, ^l^J^] another lower 
lying minimum comes up. As a typical example we con- 
sider the celebrated nucleus ^°Ne (since we take into ac- 
count one kind of particle, this nucleus corresponds to 



system with ten particles). For = the global min- 
imum occurs for a prolate shape with the configuration 
(Si,I]2,S3) — (7,7,11). With increasing the mini- 
mum initially continues to be global and the configura- 
tion remains constant while the nucleus becomes triaxial. 
However, for f2 « 4.55 another lower minimum occurs 
with a completely different, probably non-physical con- 
figuration ((Si, S2, S3) = (5, 5, 25)). Only if the original 
configuration is enforced, the oblate shape occurs at the 
larger value f^^^^^ « 5.65, but, being a local minimum, it 
may have to be interpreted as an isomeric state. 




20 30 40 50 60 70 
particle number 




70 80 90 100 110 
particle number 



FIG. 1. Mean square radii (x^), {y'^) and (z^) versus par- 
ticle number for minimal energy configurations without rota- 
tion in the triaxial harmonic oscillator. The magic (spheri- 
cal), axially symmetric and triaxial configurations are clearly 
discernible. 

The observation, that the onset of an oblate shape is 
attained at a local and not at a global minimum, pre- 
vails for all prolate nuclei with iV < 56. When looking 
at Fig.l there are the nuclei which have, for ~ 0, a 
triaxial shape close to the prolate shape, for instance 
N = 22,44,48,50. Initially when Q. is switched on, 
they behave just like the case considered in the previ- 
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ous paragraph in that they reach eventually, at ^2^^.;^, the 
oblate shape at a local minimum. For a certain value 
fi^Jj^ > the onset of tilted rotation arises. In other 
words, enforcing the initial configuration and thus per- 
taining to the local minimum, the minimum is now found 
at > and possibly (f) > while the minimum pre- 
vailing within the interval [^^crif ^critl tor 6 — has, at 

(2) 

Q. = ^^criti become unstable, in fact it turns into a max- 
imum in the 6 — (j) plane. In the following section we 
derive analytic expressions for ri^Jjt and f^^rit' P''°" 
vide an explanation why the tilted rotation cannot occur 
when starting with a prolate nucleus irrespective of the 
minimum being local or global. In this context we stress 
that the model has of course its limitations and must be 
used with due discretion. The occurrence of new global 
minima for increasing values of f2 is an expected feature 
of the model which does not necessarily have a physical 
significance. For light prolate nuclei, where the onset of 
oblate shapes occurs at values of Vt so large that the min- 
imum is not global, only a more sophisticated approach 
could ascertain as to whether such a minimum is sta- 
ble. Likewise, there is a temptation to associate the new 
global minimum with fission, but this also lies beyond 
the scope of the present paper. 

In contrast, for all oblate (and heavy prolate nuclei) 
the minima at f^^rit are global as the critical values are 
sufficiently small. In these cases we therefore view our 
findings as reliable and physically significant. As is shown 
in the following section analytically, tilted solutions can 
occur at fi^it ^ ^crit ^^"^j ™ ^^"-^ ^° occur, when starting 
from a nucleus that is triaxial for il = 0. In turn, when 
starting from an oblate nucleus, a tilted solution never 
occurs. Yet, the pattern is very similar in both cases for 



Vt < r^CTit ' ™ iiiat once the nucleus has become triaxial for 
> 0, it attains the oblate shape for Q = ri^itj where, 
for the energy minimum, the rotational axis coincides 
with the symmetry axis. 

IV. BIFURCATION POINTS 

The first subsection is devoted to the basics of our 
approach with particular emphasis on the extension of 
established knowledge. In the following subsections the 
calculation and the role of critical frequencies are dis- 
cussed. 

A. Basics 

The two critical frequencies, ^I'^^ilt ^crif closely 
related to an instability of the mean field. When H. 
is sweeping beyond the critical value, the mean field 
changes and may lead to a symmetry breaking. The 
instability can be identified as a bifurcation point pl| ] 
where new solutions that break the symmetry emerge. 
As is demonstrated below, they are related to collective 
vibrations carrying nonzero angular momentum aligned 
along the symmetry axis. The bifurcation points are de- 
fined by the vanishing vibrational frequencies in the ro- 
tating frame. That occurs at a discrete set of cranking 
frequencies given by 

fierit = UJk/K (9) 

where K is the angular momentum along the symmetry 
axis carried by a collective vibration with the energy ujk 
in the laboratory frame. 

We treat the collective excitations in the random phase 
approximation (RPA) . The variation in the one-body po- 
tential around the equilibrium shape determines the ef- 
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fective quadrupole-quadrupole interaction p^ . As a re- 
sult, the total Hamiltonian can be presented as 

2 

Hrpa = Ho-n-L-^Yl ^i^M = H-n-L (10) 

where the quadrupole operators Z?^ = r^Y2fi are ex- 
pressed in terms of the double-stretched coordinates qi = 
ZJ^lii {li — ^iVt^) (see also Appendix C). Note that 
the effective interaction restores the rotational invariance 
of the Hamiltonian in that now [H, Li] — {i = 
x,y,z). Also note that the self-consistency conditions 



Eq.(ph fix the quadrupole strength k 



5 <r^> 



T^, where 



< >=< 



y 



We solve the RPA equation of motion for generalized 



coordinates X\ and momenta Vx (see for details |23 ): 
[^fRPA,^A] = -tuJxVx, [HnPA,Vx] = iioxXx, (11) 

Here, lox is the RPA eigenfrequency in the rotating 
frame and the associated phonon operator is Ox = 
{Xx — iVx)/V^- In contrast to the CRPA approach 
P,[7[[2^ developed for the PAC rotation, the phonon is 
in the present model a superposition of different signa- 
ture phonons. The degree of the mixture depends on the 
tilted angle: the signature and \K\ are good quantum 
numbers, respectively, for rotations coinciding with one 
of the principal axes. The solution of the RPA equations 
leads to the secular equation F{lux) — with F being a 
five by five determinant. The roots are the RPA eigen- 
frequencies wa-. The non-zero solutions appear in pairs 
±hujx, we choose solutions with positive norm. The func- 
tion F factors into positive and negative signature parts 
for PAC rotations, since the signature is then a good 
quantum number. This simplifies the calculations of the 
RPA modes and provides the analytical solutions and 



the critical frequencies ^^^J^ and ^'"^l^, where the mini- 
mal solution corresponds to the PAC rotation around the 
symmetry axis being the a:-axis. 

Since at this transition point the minimal solution cor- 
responds to a one-dimensional rotation {il = (ri,0,0)) 
around the symmetry axis, the angular momentum is 
a good quantum number. The RPA states are charac- 
terised by the projection of the angular momentum upon 
the symmetry axis (x-axis) because [L^, 0\] ~ AO| with 
A being the value of the angular momentum carried by 
the phonons along the x-axis. We thus obtain 

[i^RPA,ol] = [i/-r!i.,ol] = 

iux~Xn)Ol=iuxOl. (12) 

This equation implies that at the rotational frequency 
flcr = wa/A one of the RPA frequency vanishes. 

B. Critical frequencies for oblate rotation 

1. Mean field 

Let us first turn to our numerical findings relating to 
oblate and/or near oblate cases. In these cases the first 
critical rotational frequency is associated with a one- 
dimensional rotation. We use Si < S2 < S3. Exploiting 
that at the transition point ujy — uj^ — oj^, we obtain 
from Eq.(H) a third order equation for u = Vl/uj±^ which 
reads (see also ]l8[) 



u 1 r — 1 







(13) 



2 2r + 1 

and r = S3/S2. From the discriminant of Eq.(^3|) we 
obtain the critical value 

^/27 + ^/2 

= VFWI- ^'^^ 

It was shown in [|l8| that a prolate system becomes even- 
tually oblate with the rotational axis coinciding with the 



symmetry axis, if r < rcr- In general, for r < r^r, we 
have the three solutions for Eq.jl^) 



/3r-l 

cosx = -3W- — — , 7r/2 < X < TT. 



(15) 
(16) 



2r + 1 

By their derivation these values correspond to three 
bifurcation points. We recall that for one-dimensional 







c^±/%/2 


X = 7r/2 (r = 1) ^ < 


r!2 = 









-t^±/V2 






v/273w±/2 


X = TT (r = Tcr) ^ < 




y273t^±/2 




. f^3 = 





rotation the mean field equation Eq.(13) determines the 
positive signature solutions. At the left and right hand 
boundary values of the angle x we obtain from Eqs.(p^, 

® 



(17) 
(18) 



The natural question arises: what is the meaning of the 
negative rotational frequencies? Below we demonstrate 
that the second solution gives the critical point ilc^'' at 
which the lowest positive signature vibrational frequency 
(in the rotating frame) tends to zero. 

It should be noted that the solutions 5^2,3 (u2, 3) have 
been found in ]l8[ |, however, an interpretation was not 
provided due to the lack of the RPA analysis. In j24j and 
p5[ vibrational excitations for one-dimensional rotations 
around the x-axis have been dealt with by the RPA. How- 
ever, in ||2^ the self-consistent condition and the residual 
interaction are different from those of the present paper, 
and in |^5[| the RPA analysis was done in the laboratory 
frame. We emphasise that it is important to treat the 
RPA equations in the rotating frame, since only then the 
relation between the bifurcation points of the mean field 



and the RPA solutions can be clearly worked out. This 
relation has not been analysed in the quoted papers. 

2. Positive signature solutions 

As noted above, the determinant F{ijJ\) factors into 
positive and negative signature blocks for the oblate rota- 
tion. With the projection A on the symmetry axis x of the 
angular momentum being a good quantum number 
the phonon states with different A decouple. Therefore 
the secular equation for the positive signature phonons 
separates into three independent equations each deter- 
mining the solution for the collective excitations with 
A = and A = ±2. We focus our analysis on the positive 
signature excitations with A = ±2, since they will define 
the critical rotational frequency. 

Taking into account the self-consistent condition (H) 
the secular equation can, after some tedious algebra, be 
presented for A = +2 in the form 
a be 



= 1 



2 + x 2 — X X 

using the notation 

oJx + 'lQ. 1 , r r - 1 

X— , a = r, 0= c 



(19) 



r + r r + r 



r -f 1 
(20) 



The RPA eigen-frequencies lo\ in the rotating frame are 
obtained by solving the cubic equation 



- 2a; 4c = 



leading to the solution 



(1) o /2 ^00 
uj\ , = 2lo\\ — cos 2il 

V 3 3 



(21) 



(22) 



(2,3) 

A=+2 ^ V 3 V 3 



- [cosl T %/3sin^) - 2f7. (23) 



We consider only the first solution as it is the one with 
a positive norm. Setting this frequency equal to zero 



we find just the first bifurcation point ui = Qi/luj_ of 
Eq.® 



2 X 
- cos — 

3 3 



~(i) 
A 



wliere 



~(i) o /2 X 



(24) 



(25) 



For X = 7r/2 (r = E3/S2 — 1) we obtain from 
Eq.(p5|) the well known estimate of the isoscalar giant 
quadrupole resonance energy in non-rotating axially de- 



formed nuclei, i.e. a;^"2+2 = 



with 



being 



yields 



the confinement frequency perpendicular to the symme- 
try axis of the mean field potential. The value x = tt 
2/3aj_L. For a fixed configuration 
Si < E2 < S3 the oblate state has the maximal an- 
gular momentum L^, = E3 — S2 (ground-band termina- 
tion state) HJ. If the signature symmetry is conserved 
the mode ti;^^2+2 '^^^ make a transition to a different 
configuration with two more units of angular momen- 
tum than that of the oblate (yrast) state, assuming that 
the rotational frequency exceeds the critical frequency 
rii. Therefore, for a fixed configuration the first bi- 
furcation point is related to the quadrupole vibrational 
mode w^^^2 ill ^^'^ laboratory frame carrying two units 
(A = +2) of angular momentum. 

For the A = — 2 excitations the secular equation is 
similar to that of Eq.(p^ with x replaced by —y where 



y 



UJ±_ 



(26) 



The RPA eigen-frequencies lu\ are obtained by solving 
the cubic equation for y 



y 



2?/ - 4c = 



(27) 



i"'^^2 = 2il + cj^ A fcos ^ ± VSsin ^ 



(28) 



4!l_2 = 2r!-2c.iy^cos|. (29) 



The positive norm solutions are the first and the third. 

(1 31 

Setting uj^^_2 equal to zero, we obtain the second and 
third bifurcation points U2,3 — ri2,3/^± of Eq.(p^) which 
read 



/2 

2 V 3 V 3 



cos ^ + 



~(3) 

A 

1) 

A=-2 



~(i,3) /2 / X . . X 

uj\_ \ ^LO^d - [cos — ± v3sm — 



(30) 
(31) 
(32) 



~ (1 3) 

The frequencies uj\^_2 are the vibrational modes in the 
laboratory frame carrying A = — 2 units of angular mo- 
mentum. 

A collective vibration with the projection oi K ^ Q 
units of angular momentum along the symmetry axis may 
be associated with a surface wave travelling around this 
axis. This wave can travel in opposite direction (negative 
rotational frequency) with respect to the rotation of the 
core (the vacuum) 

With X = 7I-/2 (E3 = S2) we obtain CjxI-2 = V2t^± = 
'^i=+2 ^'i'^ ^'x=~2 ^ ^- ^^i^ case the angular mo- 
mentum of the vacuum (oblate) state is zero, since — 
S3 — S2 = 0. Thus, in the oblate system with zero angu- 
lar momentum a breaking of the symmetry arises at the 
positive/negative rotational frequency 



(33) 



It is at this point where the transition from non-collective 
rotation (around a;-axis) to triaxial collective rotation 
takes place, where the latter state has two units of angu- 
lar momentum. 



leading this time to the solutions 



At the other end of the parameter range, x = tt (r 



Tcr), we obtain cj^"2-2 ~ 2-y/2/3( 



and a; 2 



— •\/2/ StJ j_ • Within the fuU range of parameter values 
7r/2<x<7r(l<r< rcr) the vibrational modes uj^^_2 
are sweeping over the values 

^ (34) 
(35) 



0. 



When the rotational frequency approaches the lowest 
critical value for a triaxial nucleus which is 



no 



2 / X 
(cos- 



V3sin|) >0, (36) 



the vibrational excitation in the rotating frame for the 
oblate state tends to zero, i.e., w^t^ j = '^i=-2 + 



2^21 



~,(3) 



2n 



(1) 



0. The negative 



quadrupolc mode energy wj^^ ^ vanishes at the transi- 



tion from the non-collective to the collective triaxial ro- 
tation and becomes positive for the oblate state at higher 
rotational frequencies of the mean field. In other words, 
the mode cjj^^ ^ can be viewed as a de-excitation in the 
intrinsic frame of the vacuum (oblate) state thus leading 
to the triaxial state with two units of angular momentum 



(1) 

A=-2 



less than the vacuum state. In contrast, the mode uj 
is positive at S^^rit ^^"^ describes an ordinary quadrupole 
excitation with two units of angular momentum less than 
that of the yrast state on which it is based. 



3. Negative signature solutions 

From the previous subsection we have learned that 
the critical rotational frequency is related to the van- 
ishing of the RPA solutions in the intrinsic frame and 
consequently to the breaking of the mean field symme- 
try. Since the negative signature solutions are associated 



with the A = ±1 excitations, at least one of the solutions 
must be related to signature symmetry breaking, i.e. to 
the bifurcation point where tilted rotation is setting in. 

In the oblate case the secular equation for the negative 
signature excitations contains the "spurious" solution at 
Lu — Q corresponding to a collective rotation and the vi- 
brational modes; the latter include the wobbling motion 
and carry one unit of angular momentum [^,0. Taking 
into account the self-consistent condition (||) the secular 
equation for the A = — 1 excitation reads 

/ 



- uj± + — X ijJx + '-^^± — Q + X 
d p 



UJx — LU± — Q + X UJx — LU± + Q — X 

where we use the notation 

1 + fc „ r + k 



2 



(37) 



r+ 1' 



/ 



a = , p = 



r + r 

r — k 
r + 1' 



k = E1/S2. 



The root in x of Eq.(|37[) is the RPA frequency uj\. After 
some simplification we obtain 



y{y^ -V- 2ac) = 0, 



where 



y 



Ll!\ — Q 

V^x + ^: 



(38) 



(39) 



The trivial solution y — Q corresponds to the "spuri- 
ous" excitation lj\ — Q in the laboratory frame with the 
mode LJx = u}\ + Q being a rotational mode with fre- 
quency VL in the rotating frame. For Eq.(^8|) there exist 
three other solutions 

ij} -f 27rn 



uJn = 2\ — - cos ■ 



cos t/i 



27- 



3 

r 



1 



n = 0,1,2 (40) 
, 0<^<7r/2 (41) 



(cc;2+cj2)3/2^ + l^ 
The first solution (n = 0) is the one with the positive 
norm. Setting it equal to zero, we obtain the first bifur- 
cation point for the negative signature 
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^^.=-2/5+^cos| = % (42) 

Here, the mode 

^iU-2^^^oost (43) 

describes the quadrupole excitation in the laboratory 
frame, which is based on the yrast (oblate) state and car- 
ries one unit of angular momentum, i.e. A = — 1. Since 
cos J > \/3/2, fli is larger than the maximal rotational 
frequency (see Appendix B), it cannot be associated with 
a stable tilted rotation. For S3 = S2 (V' = 7r/2) the mode 
oj^xl-i ^ V ^ again the isoscalar quadrupole gi- 

ant resonance energy in non-rotating axially deformed 
nuclei |^^. Notice that this mode has been introduced 
in the context of a BEC under the name "scissor mode" 
|p7| . In the spherical case it is w^^-i = V^(^o ||]- With 
V' = it is = 2y/ {luI+ujI)/3. 

For the A = 1 excitations the secular equation is similar 
to Eq.(|37|) with x replaced by —x. After simplification, 
we obtain 

yiy^-y + 2ac)^0, y = ^^±£^. (44) 

Again, the trivial solution y = corresponds to the " spu- 
rious" excitation a)> = in the laboratory frame with 
the mode lj\ = ujx — VL being a rotational mode with 
frequency —51 in the rotating frame. The three other 
solutions of Eq. ( ^4| ) are 

- (cos I T V3 sin I) ™ (45) 

-E,. = -2y^^cos|-n. (46) 

The positive norm solutions are the first and the sec- 
ond. Setting equal to zero, we obtain the bifurca- 

tion points 



.../5^(cos|-V3sm|).% (4, 

^^3^/^^3(cos| + V3sin|).i% (48) 

4t;l-/^^(cos|TV3sm|). (49) 

~ (1 2) 

Thus, they are related to the vibrational modes in 
the laboratory frame and carry A = units of angular 
momentum. 

In fact, the lowest mode w^^^j^(il) corresponds to a 
wobbling excitation in the intrinsic frame. As soon as 
the rotational frequency approaches the lowest value ri2, 
the mode ~ '^a=+i ~ ^ becomes soft and signals 

the onset of genuine tilted rotation of the nucleus. The 
mode energy vanishes for 

^. = 2/^^3cosi±^ = f^^?), (50) 

Thus, the bifurcation frequency U,2 corresponds to the 
transition from dynamical fiuctuations of the angular mo- 
mentum in the oblate system to the stable tilted rotation 
of the triaxial shape. This critical frequency coincides 
perfectly with our numerical results. 




1 1.05 1.1 1.15 1.2 1.25 1.3 
£2 



FIG. 2. Critical frequency il'^'^ versus r — E3/E2 for a 
typical (near) oblate situation. The example refers to A'' = 16 
with (El, E2, E3) = (30, 40, 44); with lj^, = 20.7, uji_ = 14.5 it 
is VL^^X = 0.93 for r = 1.1. 

Furthermore, Eq.(pO|) explains why oblate nuclei (at 
= 0) cannot make the transition to tilted rotation. 
In fact, for these nuclei r = 1 and hence i]j = 7r/2 in 
Eq. (^ . As a consequence, f2^^;\ is zero (ci)^2+i = 0)i 



i.e. no transition can occur with our convention Si < 

E2 — S3. In contrast, whenever r > 1, a positive vahie 

(2) 

is obtained for il^/t as illustrated in Fig. 2. 

With r = 1, i.e. L-c = S3 — S2 = 0, the other two modes 

~,(2) 



A = ±1 are degenerate, since (^x=+i = V^x + 
io^^l-i- T^^^ corresponding critical frequencies 



(51) 



are larger than the maximal frequency (see Appendix 
B) and do not give rise to symmetry breaking of the 
mean field. Nevertheless, these critical frequencies can 
be related to vibrational excitations as discussed in the 
following. 

As was suggested by Marshalek and Sabato the 
self-consistent cranking model should allow the study of 
vibrational states in spherical nuclei rotating around the 
symmetry axis. Our results for quadrupole excitations 
demonstrate that this idea can be succesfuUy employed 
for any axially deformed system rotating around the sym- 
metry axis (see also a general discussion in pl| ). The 
vibrational excitations with nonzero angular momentum 
A (or K) aligned along the symmetry axis are determined 
by the bifurcation points of the mean field via Eq.(|^) 



0J\ — Ailcr- 



(52) 



The lowest critical rotational frequency (bifurcation 
point) gives the collective excitation with the largest an- 
gular momentum along the symmetry axis. If the crit- 
ical rotational frequency lies below the maximal rota- 
tional frequency, it will give rise to a mean field symme- 
try breaking, as the corresponding vibrational excitation 
vanishes in the rotating frame. If the frequencies lie above 
the maximal value, they yield the higher vibrational ex- 
citation energies with nonzero angular momentum in the 



laboratory frame. They persist as their counterparts in 
the rotating frame do not vanish. 

We have established numerically in all cases of triax- 
ial oblate nuclei the onset of tilted rotation and iden- 
tified exactly the value of Eq.(|50|) as the lowest rota- 
tional frequency. The values of are all fairly small 
in these cases (f^^u < 1)' since the values for r are close 
to unity. This is in contrast to light prolate nuclei to 
which we turn next. For prolate nuclei (at = 0) we 
have Si = S2 7^ S3 and the oscillator frequencies are 
u)^ — ujy — u)'^ implying the symmetry axis to be the 
z-axis. Since in this case r > 1, Fig. 2 cannot serve as a 
basis to exclude tilted rotation. Yet we present an inde- 
pendent argument. 

C. Critical frequencies for prolate rotation 

Consider a rotation around the symmetry axis, i.e. the 
z-axis. The eigenmodes have a particular simple form 
p"5|-p7[ and Eq. (|^) leads to the nontrivial solution which 
must satisfy the equation 



(w+cj_ - n'^) = 0. 



(53) 



Setting ujx = LOy = uj^ we obtain the bifurcation point 



(54) 



For > r^cr the axial symmetry is broken and the sys- 
tem is driven into the domain of triaxial shape under the 
PAC rotation. This situation is reminiscent of a striking 
feature established experimentally in the rotating BEC 
32[ ; there the minimal rotational frequency for nu- 
cleation of a vortex occurs at around 0.7luj_, where luj_ 
is the mean oscillation frequency of atoms in the x — y 
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plane, irrespective of the number of atoms or oscillation 
frequency along the z-axis. According to our analysis, 
the dynamical instability in a nuclear system and in the 
rotating BEC is of a similar nature, in spite of the dif- 
ferent character of the interaction between nuclcons (at- 
tractive) and between atoms (repulsive). The important 
aspect is the trapping of either system by the harmonic 
oscillator potential. For the nuclear system it is the mean 
field and for the BEC it is the external magnetic field. 

For a prolate system the projection K of angular mo- 
mentum is a good quantum number as [Lz, Oj^-] = KO^^^. 
We thus obtain an equation similar to Eq.(|l2|) after re- 
placing X hy K . At the rotational frequency ficrit = 
Cok/K one of the RPA frequency vanishes. For a pro- 
late rotation the RPA equations are solved separately for 
each quadrupole mode, since there is no coupling between 
different K modes. 

For the quadrupole phonons with the largest projection 
if = 2 it is 

uJK=2 = V2uj^ - 2VL. (55) 

This mode is the quadrupole excitation having two more 
units of angular momentum than the vacuum state. 
When ujk=2 — the transition from non-collective ro- 
tation (around the z-axis) to triaxial collective rotation 
takes place, i.e. at 

f]„ = ^ (56) 

which is just the bifurcation point of the mean field in 
Eq. (^4[) . This bifurcation point applies for any axially 
symmetric system (prolate and oblate) as follows from 
Eqs.(|3|) and (||). 

For the K ^ \ mode the RPA solution is 

ujK=i^ J^^+^l-^- (57) 



The condition ijJK=i = yields the critical frequency 
at which the onset of the tilted rotation should occur for 
the prolate system, i.e. at 

n„ = yC^T^. (58) 

However, the energy to create this mode of angular 
momentum l?i is too high, the system rather prefers 
the PAC rotation around the axis perpendicular to the 
symmetry axis. Besides, as discused above, the value 
of the rotational frequency in Eq.(p8|) lies outside the 
physical range of rotational frequencies (see Appendix 
B). As a consequence, our model does not allow tilted 
rotations for systems which are prolate at = 0, i.e. if 
Til — "2,2- In contrast to the rotational frequency defined 
by Eqs.(^), (p6|), this frequency is related to vibrational 
excitations carrying one unit of angular momentum in 
an axially symmetric system, be it oblate (Eq.(|5l|)) or 
prolate (Eq.(|5|)). 

Note that the assumption Si = E2 is essential for the 
arguments presented here. In contrast, if all three val- 
ues of Yii differ from each other, tilted rotations do occur 
at the value given by Eq.(^o|) just as it has been found 
for the near oblate situations. We recall, however, that 
this happens for < 56 at values so high that the cor- 
responding energy minimum is only local. In contrast, 

(2) 

when A^ > 70, the relevant value of is sufficiently 
small so as yield a global minimum. 

V. SUMMARY 

Within the model considered only nuclei being triaxial 
in their ground state exhibit tilted rotation at a specific 
rotational frequency. In fact, the analytic expressions for 
these critical frequencies obtained in the present paper 



are in a perfect agreement with our numerical findings. 
All nuclei attain at ^i^^^^ an oblate shape. The value of 
^^^11^ depends on the particular nucleus. For light and 
medium prolate or near prolate nuclei {N < 56) this 
energy minimum is a local minimum, while for oblate or 
near oblate nuclei and heavy (near) prolate {N > 70) 
it is a global minimum. All triaxial nuclei - these are 
the nuclei which form the little 'bubbles' in Fig.l ~ go 
through a second transition and admit a tilted rotation 
beyond f^^rit "while the tilted rotation cannot occur when 
starting with a nucleus that is axially symmetric at Q — 
0. From our model, the occurrence and identification of a 
tilted rotation must be seen as a signature for triaxiality. 
It is true that the harmonic oscillator model is a simplistic 
approximation for the phenomenon of the tilted rotation. 
However, our findings may provide guidance for values of 
particle number and rotational frequency at which tilted 
rotations may be found. 

Using the cranking approach to describe non-collective 
rotation around a symmetry axis we have established the 
correspondence between symmetry-breaking bifurcation 
points of the mean field and the RPA frequencies carrying 
non- vanishing angular momentum aligned along the sym- 
metry axis. We have demonstrated that the dynamical 
instability in a nuclear system and in the rotating BEC 
are of a similar nature, in spite of the different charac- 
ters of the interaction between nucleons (attractive) and 
between atoms (repulsive). The remarkable agreement 
with regard to the bifurcation points in both systems is 
due to the common nature of the confining potential and 
its dependence on the rotational speed. 
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APPENDIX A: EIGENMODES 

The Hamilton function of Eq. (|l|) can be written in ma- 
trix form 



H = {p,r\^n{p,r^ 



(A.l) 



where p and f are combined to the six dimensional col- 
umn vector {p, f} and 



n = 



< I 








1 












1 










-^z 






-a 





2 



(A.2) 



We aim at the quantum mechanical form in terms of bo- 
son operators 

H = {Q,Q^^n,^{Q,Q^ = Y,E,{Q]Q, + 1/2) 

(A.3) 

where we denote by {Q, Q^} a column vector which is 
the transpose of the vector (Qi, Q2, Qs, Q3, Q21 Qi) 
where 



TLn 



( 






E2 
\Ei 















^2 

5:3 









Ei\ 





/ 



(A.4) 



It is instructive to obtain the eigenmode energies from 
a purely classical calculation. This approach also pro- 
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vides the linear transformation between the Qj and the 

{Pt,Xi}. 

The classical equations of motion for the Cartesian 
components of the momentum and position coordinates 
read 



dt \r 

The matrix M is given by 



M 



(A.5) 



(A.6) 



where /a is a 3 by 3 unit matrix. The classical orbits are 
the solution of Eq. ( [A.5| ) which reads 

(A.7) 



r(i)} = U exp(I?t)V{p(0), f(0)} 



where V = dia.g{—iEi,—iE2,—iE3,iE3,iE2,iEi) is the 
diagonal form of M. The initial conditions {p(0),r(0)} 



are of no interest here. From Eq.(AJ) it follows that 



the eigenmodes listed in Eq.(A.4) occur in the classical 
problem in this form. They are obtained from the secular 
equation 



det|£;/6 = 0. 



(A., 



The above equation turns out to be a third order polyno- 
mial in E'^ and has also been found by IHl. The column 
vectors of U are the (complex) right hand eigenvectors 
of M. Note that M is not symmetric, hence neither U 
nor V are unitary, yet V —lA^^. We denote the column 
vectors of U by u'^^\ they obey the equations 



(A^ + iSfc/)!*^^) = 0, /c= 1,2,3 
{M-iE^_kI)u'^^^ ^Q, A: = 4, 5, 6 



(A.9) 



which can be solved as an inhomogeneous system by 
choosing an arbitrary component of u'^^^ equal to unity. 
The proper normalisation is achieved by the observation 



that, up to normalisation factors, the matrix V = Li ^ 
can be written as 



V 































—i 














—i 














i 














i 




























0/ 





h 



-h 





(A.IO) 



This implies that U^HU is in fact skew-diagonal as in 



Eq.(A.4), and therefore U can be normalised such that 
14^1-114 = Ti.qm ■ This implies in particular 



{p,r}^U{Q,Q^}. 



(A.ll) 



The properties of the matrix U guarantee the com- 
mutator [Qk,Qki] = Sk,k' as a consequence of [xi,pj] = 
ihSij. Clearly, for = the Qj become the usual oper- 



ators tti — (j>i — imujiXi)/ ^2mhuji. 



APPENDIX B: MAXIMAL ROTATIONAL 
FREQUENCY 



At particular values of the roots of cubic equation 
become negative, i.e., the corresponding eigenmodes E^. 
would be imaginary. Let us consider the simple case of 
a rotation about the principal axis x. The solutions for 
the eigenmodes are well known p5|-p^ 



El — E2 — -E3 — ^~ 



„ loI+ujI „ 1 
jI = y ^ ^ + ± - 



(B.l) 

1/2 



(B.2) 

Thus, for the eigenmodes to be real, ujy, lo^ and must 
satisfy 

-1 1/2 



(B.3) 
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which leads to the inequahty 



(B.4) 



merits using x-axis quantization, since the symmetry axis 
coincides with the axis of rotation. 



This means that the frequency must be either larger 
or smaller than both frequencies. Since the condition 
U, > sup(ci;j,, ujz) is unphysical, the rotational frequency 
should not exceed the maximal rotational frequency 

f^max < inf(Wy, UJz)- 

APPENDIX C: QUADRUPOLE OPERATORS 

For a solution of the RPA equation of motion it is 
convenient to use the following representation of the 
quadrupole operators 



Do 
Di 

D2 

D3 



I6n 



{2f 



f) 



V2 



{D2 



D2-1) 
D2-2) 
D2-1) 



-yz 



-j={D22-D2-2) 



15 ,_2 _2x 

i6^(" -y^ 



(C.l) 
(C.2) 
(C.3) 
(C.4) 
(C.5) 



With these definitions we have Y^n.=^2 -^u-^m ^ 



Sm=o^ni- Using the transformation Eq.( A.ll ) the 
quadrupole operators are expressed in terms of the boson 
operators Qi,Q\. In case of the one-dimensional rotation 
the matrix hi reduces to the transformation as in and 
all matrix elements are calculated analytically. When the 
signature is a good quantum number, the operators can 
be classified according to their signature quantum num- 
ber: the operators Dq, Di, D2 are characterized by the 
positive, while the operators D3, D4 are characterized by 
the negative signature quantum number |^,|^,^. In the 
oblate case it is convenient to calculate the matrix ele- 
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